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Fast amplitude variations due to atmospheric scintillation are the main concerns 
for the Deep Space Network (DSN) Ka-band downlink under clear weather con- 
ditions. A theoretical study of the amplitude scintillation variances for a Unite 
aperture antenna is presented. Amplitude variances for weak scattering scenarios 
are examined using turbulence theory to describe atmospheric irregularities. We 
first apply the Kolmogorov turbulent spectrum to a point receiver for three dif- 
ferent turbulent profile models, especially for an exponential model varying with 
altitude. These analytic solutions then are extended to a receiver with a finite 
aperture antenna for the three profile models. Smoothing effects of antenna aper- 
ture are expressed by gain factors. A group of scaling factor relations is derived to 
show the dependences of amplitude variances on signal wavelength, antenna size, 
and elevation angle. Finally, we use these analytic solutions to estimate the scin- 
tillation intensity for a DSN Goldstone 34-m receiving station. We find that the 
(rms) amplitude fluctuation is 0.13 dB at 20-deg elevation angle for an exponen- 
tial model, while the fluctuation is 0.05 dB at 90 deg. These results will aid us in 
telecommunication system design and signal-fading prediction. They also provide 
a theoretical basis for further comparison with other measurements at Ka-band. 


I. Introduction 

Turbulent irregularities exist in the Earth’s lower atmosphere — the troposphere. These mixed air 
masses act as moving parcels with various sizes, randomly varying in space and time. The parcels of 
air usually have different temperature, density, humidity, and pressure across their spans. Thus, there 
are different refractive indexes across them. When radio waves pass through these turbulent cells, their 
propagation directions and path lengths are altered due to diffraction and scattering. Rapid variations 
of signals in phase and amplitude thus occur. The wavefront becomes tilted and crinkled in a plane 
perpendicular to the propagation path. As a result, received signals from all directions and for all 
instantaneous components are superimposed. This results in the signal being enhanced or diminished 
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in both phase and amplitude. Random fast fluctuations around the average of field strength are called 
amplitude scintillations and are the subject of this article. 

It is generally known that scintillation intensity (also called fading depth) and signal amplitude increase 
significantly with increasing signal frequency. This is because the atmosphere has many more turbulent 
cells with small sizes than larger ones. When the wavelengths of radio signals match the Fresnel lengths 
of the cells, the scintillation effects become intense. Thus, 32-GHz (Ka-band) signals are expected to 
have larger fluctuations than those around 8.4 GHz (X-band) and around 2.0 GHz (S-band). The fading 
depths also increase rapidly with decreasing elevation angles. Because a low elevation path corresponds 
to a longer atmospheric column, more turbulent cells are included in the propagation path. However, the 
antenna size tends to reduce scintillation effects since the influence of small turbulent cells is averaged 
out by large apertures. 

Atmospheric scintillations have been extensively studied in both experiments and theories through the 
last half century [1-4]. A variety of scintillation models have been developed. Turbulence theory becomes 
the essential in understanding refractive irregularities in the lower atmosphere. The Kolmogorov model is 
able to explain the most important features of turbulence which cause the perturbation of radio waves. In 
order to increase the transmission bandwidth, NASA is going to upgrade its Deep Space Network (DSN) 
downlink frequency to Ka-band. Under fair weather conditions, amplitude variations due to atmospheric 
scintillation will become the main concern most of the time. This is because fast and deep fadings can 
cause receiver loss of lock. Under clear weather, at Ka-band the dominant perturbation on receiving 
signals will come from the scintillation effects. We need to understand how severely this will impact on 
amplitude variances and how fast the fading can be caused at Ka-band. 

In this article, we will revisit turbulence theory to understand how the atmospheric turbulences cause 
the scintillation for an antenna with finite aperture. Such a detailed investigation for a large antenna at 
Ka-band has not been done, primarily due to the lack of experiments and applications. For example, if 
a large antenna aperture can smooth the influence of small turbulent eddies, how much will amplitude 
variances be reduced? We will start from the random wave equation to develop an analytic solution for a 
point receiver. Then we will apply the solution to three profile models to describe the height-dependent 
variation of atmospheric turbulence and then integrate each over a large antenna area. We will apply the 
Kolmogorov turbulent spectrum into a DSN Ka-band downlink scenario to calculate amplitude variances 
for the Goldstone receiver. We will re-examine the dependences of amplitude scintillations on frequency, 
elevation angle, and antenna size for DSN receivers. This study will provide a theoretical basis for any 
further comparison study with Ka-band measurements. We also extend this study to an exponential 
profile model, which has not previously been explored. 


II. Turbulence Theory 

When radio waves pass through the turbulent atmospheric media, the instability of the electric field 
amplitude is affected by the medium’s dielectric constant variance (assuming that the atmosphere has 
a homogeneous composition, but inhomogeneous temperature and pressure). The relation can be found 
through solving the following random wave equation: 


V 2 E+k 2 [l + 6e(r,t)\E = 0 (1) 

This results from Maxwell’s equations for the electromagnetic field [1, Section 2.1, pp. 6-10]. Using the 
R.ytov approximation to solve this equation for a small Se perturbation in the dielectric constant e, the 
electric field strength expression becomes [2, Section 2.2, p. 14-19] 
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(2) 


E(R) 


Eq(R ) exp 



J d 3 rS£(r)A(R,r) 


where 


A(R,r) = 3? 


G(R,r) 


Eo(r) 

Eq(R) 


( 3 ) 


is the real part of the Green function G(R,r ) and the ratio of electric fields. We emphasize the natural 
logarithmic of the amplitude fluctuation: 


X = In 


E 

Eq 


Fluctuations of the field strength A E are small compared to Eq since microwave propagation is described 
as weak scattering: 


Eq + AE\ 

Eq J 


= In 



A E 
~Eq 


< 1 


( 4 ) 


Thus, from Eq. (2), 


X = ~k 


2 


d 3 r5s(i')A(R, r) < 1 (Np) 


( 5 ) 


If converting \ into decibel units, we have utilized the following relationships: 

X'(dB) = (10 log 10 e)x(Np) = 4.34x(Np) (6) 

<X 2 ) (dB) 2 = (101og 10 e) 2 <x 2 ) (Np) 2 = 18.84 ( X 2 ) (Np) 2 (7) 

Because the amplitude fluctuations are referred to their mean value, the ensemble average (represented 
by an angle bracket) of the amplitude fluctuation vanishes: 

(X) = 0 (8) 

Its mean-square value is not zero, and we can use Eq. (5) to describe the scintillation intensity: 

(x 2 ) = k 4 J d 3 r J d 3 r' A(R,r)A l (R,r') (fe(r, t)fc(r / ,t)) (Np) 2 (9) 

Because e = n, we can write Se = 2no6n. When uq — 1, this becomes 5e = 2 Sn, and 

( Ss(r , t)6e(r', t)) = 4 ( Sn(r , t)Sn(r', t)) (10) 
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We know that, in the r < — > k domain [1, Section 2.2.5, pp. 21-23], a three-dimensional Fourier 
transform connects the spatial covariance of refractive index fluctuations B n {r ) and the wavenumber 
spectrum of irregularities $„(«;): 


B n { r,r') 


(i Sn(r , t)6n(r',i)) 


j d 3 K,$ n (n) exp [in ■ (r — r')] 


( 11 ) 


The refractive-index structure function is defined as 

V„{p) = ( [<5n(r, t) - 5n(r + p, t)]~^ 


(12) 


and is often approximated by the following approximation due to Kolmogorov: 

Vn{p) = C 2 n | P | 2/3 


(13) 


where C/ is the refractive-index structure constant and p is the spatial difference between two points r and 
r'. Previous studies usually have assumed a uniform atmosphere so that Cl is a constant independent of 
altitude. In this article, we will consider three different altitude profiles of Cl so as to better characterize 
the real atmosphere. Based on the 2/3 power law on the spatial separation, the wavenumber spectrum 
of irregularities, $„(«;) for the Kolmogorov refractivity model is expressed as 

$„(«;, z) = O.O33C^(0)k“ (11/3) , 0 < k < oo (14) 


III. Amplitude Variance for a Point Receiver 

Consider the vertically downward propagating microwave shown in Fig. 1. After encountering a turbu- 
lent eddy, the signal will be diffracted with a small angle 0. In this article, we consider only weak forward 
scattering. This means that 9 = A/l <C 1, where A is the wavelength of the radio signal (~0.01 m) and l is 
the size of the turbulent cell. Notice that the first Fresnel length is given by _//, = \f~Az. The wavelength 
of the signal is much less than the size of the turbulent eddies. Because 9 <C 1 (about 1 mrad), and r <C z 
in Fig. 1, the downcoming wave stays quite close to the vertical axis. This condition is called the paraxial 
approximation. The coordinates used in this article are shown in Fig. 2. 

For an overhead source, we use the cylindrical coordinates identified in the figure with the point 
receiver placed at the origin, R = 0: 


E(r ) = E 0 exp(—ikz) 


(15) 


E(R) = E 0 


(16) 


G(R,r) 


1 

4irp 


exp (ikp) 


(17) 


where p = \/ z 1 + r 2 . Both displacement vector ( p ) and turbulent wavenumber vector ( k ) in cylindrical 
coordinates are expressed as 
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Fig. 1. Geometry showing how a micro- 
wave passing through a turbulent refrac- 
tivity eddy is scattered to the receiver. 



Fig. 2. Coordinate systems describing the displacement vector 
p and the turbulent wavenumber vector k. 
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p = zi z + r cos (j)i x + r sin <jA y 


(18) 


k = K z i z + K r cos u;i x + K r sin u>i y 


(19) 


Thus, for the amplitude weighting function introduced in Eq. (5), we have 


A{R,r) = 3 ft 


exp (ik\J z 1 + r 2 ) E 0 exp {—ikz) 
47 tV z 2 + r 2 E 0 


(20) 


In order to get an analytic solution and in view of the small-angle scattering indicated by Fig. 1, we can 
expand the square-root distance as follows: 


, — / 7^2 r 2 

^ + ,■2 = ^1 + -^+- 

In the numerator of Eq. (20), we must keep both terms but need only the first in the denominator. The 
amplitude weighting function thus becomes 


A(R,r ) = 3? 


exp ikz + ik — 
2z 

47 TZ 


exp (—ikz) 


or 


A(R,r) 


1 

- — cos 

47 TZ 



(21) 


With this simplified expression, from Eqs. (9) and (10), we can write the variance of logarithmic 
amplitude fluctuations as 


<x 2 > 


fc 4 

47T 2 



ridri cos 




x 



r 2 dr 2 cos 



(6n(p 1 )Sn(p 2 )} 


(22) 


In spherical wavenumber coordinates, the spatial covariance of refractive index fluctuations can be ex- 
pressed as follows: 


(Sn(p 1 )Sn(p 2 )) 


d 3 n$ n 


Zi+ z 2 \ 


exp [in ■ \p 1 


Pal] 


2 , r j, • ( z l + z 2\ f , 

J k dn J dipsmip<s n Ik, — - — I J aw 

x exp {i[z 1 KCOS'i/> + r 1 Ksin^cos(</>i — cu)] — i[ziKCOStl> + riKCOsipcos(4>2 — w)] } (23) 
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where we have allowed the turbulent spectrum to vary with wavenumber and the average height of the 
irregularities. Thus, 


, 2V k 4 r dz x r°° dz 2 f°° , r°° , r 2j r,, ■ / 

<X ) = y — n / / / TWi / r 2 dr 2 / ndn dipsmip$ n h 

4tt- J 0 zi 2 2 Jo Jo Jo Jo V 


2l + 22 


f kr 2 \ ( hr 2 \ / >2lr 

x cos ( — ! - ) exp(*ZiK:cos cos ( — — ) exp(— iz 2 «:cos ip) / du> 

\22i J \zZ2 J Jo 

p2it p 27 T 

x / exp [riKsin J>cos(J>i — / exp [r 2 K; sin J>cos(J> 2 — w)]d</> 2 

Jo Jo 


(24) 


Using the integral relation Eq. (A-l) in Appendix A to integrate difiidfodui , we can reduce this as follows: 


POO P7t 

(x 2 ) = 27 tA: 4 / K 2 dK / dip simp 
Jo Jo 

f°°dz 1 dz 2 r . / , /lA ( 

x / / exp i zi -z 2 kcosw k, 

Jo 21 Jo 2 2 V 


21 + 22 


r 00 f kr 2 \ f°° f kr 2 \ 

x / ridricos ( — - I Jo(riKsin^)) x / r 2 dr 2 cos ( — - I J 0 (r 2 «:sin'0) (25) 

Jo V 22i J Jo V 222 J 

Using the relation Eq. (A-2) in Appendix A, we can perform the radial integrations to find 


P OO /*7T /»00 /»00 

(y 2 ) = 2irk 2 / k 2 <Jk / dip sin ip / cfoi / dz 2 exp [mcos ip(z x — Z 2 )] 

Jo Jo Jo Jo 


2/c 


2 i«rsin r/A f Z 2 K 2 sin ip \ ( Z\ + z 2 

x sm I — sin — I <1>„ I k 


2fc 


(26) 


The variation of the exponential term is very rapid and destructive unless the values of z x and z 2 are 
very close to each other. This lets us approximate the double integral as follows: 


pOO p OO 

/ dz\ / dz 2 exp \in cos ip(z x — Z 2 )]F (z x , Z 2 ) — 2ir5 (k cos ip) F (z, z) 
Jo Jo 


(27) 


where 

2l + 2 2 
^ ~ 2 

where we exploit the Delta function relations provided in Eqs. (A-5) through (A-7) in Appendix A to 
write 
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(X 2 )=4A- 2 J ndn J dipsmip$ n (n, z) J S (cos ip) sin 2 dz (28) 

Since the delta function requires that cos ip = 0 and sin ip = 1, we have finally 

r°° ( zk 2 \ 

(y 2 ) = 47T 2 fc 2 J KdK<$> n (K,z) J sin \~ 2 kjdz (29) 

IV. Results for Three Turbulent Profiles 

We will apply the Kolmogorov spectrum to describe the turbulent irregularities for a point receiver. In 
doing so, we consider three profile models to describe the altitude variation of turbulence strength, as rep- 
resented in the parameter C 2 . The profile models are (1) uniform slab, (2) thin layer, and (3) exponential. 
They are illustrated in Fig. 3. 

A. Slab Model 

In this case, we assume that all atmospheric turbulent eddies are of uniform strength and contained 
within a limited height, H : 


SLAB LAYER THIN LAYER EXPONENTIAL 

MODEL MODEL LAYER MODEL 

Fig. 3. Three profiles used to describe the level of atmospheric 
turbulence: (1) slab with uniform density, (2) thin layer with thick- 
ness A H , and (3) exponential decayed atmosphere. 
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After applying the Kolmogorov spectrum and slab model into the expression of amplitude 
Eq. (29), we have 


<X 2 > = 47T 2 fc 2 



r H 

ndn / dz$ n (K, z) sin 

J 0 



dz 


poo 

= 2ir 2 Hk 2 j KdK 0.033 C 2 ^ 11 ^ 


sin (Hn 2 /k) \ 

Hu 2 /k J 


= 0.651 k 2 HC 2 



dn 
K 873 



sin (Hk 2 /k) \ 

Hn 2 /k J 


Substituting £ = Hk 2 /k, we have the following for the point receiver: 


<X 2 > = 0.326Cjff"/« t VejT_^_ 



sin£ 

T" 


or finally 


(X 2 ) = 0.307 C 2 i7 11/6 fc 7/6 (Np) 2 


B. Thin Layer 

A thin layer of thickness AH often is used to describe a turbulent layer or a cloud layer 
altitude H. This is defined as [4-6] 



H < z < H + AH 
otherwise 


When AH << H, we have 


<X 2 > = ^ 2 k 2 


ndn<& n {n, z) 




dz 


pOO pC 

27 r 2 k 2 / ndn<& n {n, z) 

Jo Jo 


' / zk 2 

1 - cos — — 
V 2 k 


dz 


= 0.651k 2 AHC 2 


dn 


to «' 


8/3 


1 — COS 


Hk< 


Using a substitution of ( = Hti 2 /k, we have 

(x 2 ) = 0.326C 2 AHH 5 / e k 7/6 [ [1 - cos C] = 0.563C 2 AHH 5 / 6 k 7/6 (Np) 2 

Jo C ' 


variances, 


(31) 


(32) 


(33) 
at a fixed 

(34) 


(35) 


(36) 
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C. Exponential Model 

In the third case, we assume that atmospheric turbulence level has its maximum value at the Earth’s 
surface and decreases exponentially with altitude. It is defined as 

CnW = C'iU 0 ) ex P (-Jp) ( 37 ) 


where H is the scale height. Thus, 


<E 


t (K, z) = 0.033k (11/3) U 2 (0) exp 


We have 


( X 2 ) = 4 ^ fc 2 J°° KdKOm3 k -(H/3 )( o 2 (0) ^°° exp gin ^^0 


OO J tym pOO 


= 0.651fc 2 U 2 (0) 


dhi 


exp (-s) 


i ( ZK 2 

‘-“a 




= 0.651fc 2 C 2 (0)if [ 

Jo 
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Hu' 


« 8/3 /Hk 2n 2 


1 + 


Using a substitution of £ = Hn 2 /k, we have 


<X 2 ) = 0.326fc 7 / 6 C 2 (0)if 11 / 6 ^°° c f 1 C / 6 i f C 2 = 0.530fc 7 / 6 C 2 (0)Ef n / 6 


(38) 


(39) 


(40) 


Using the integral result, Eq. (A-10), we have 

( X 2 ) = 0.530fc 7/6 C' 2 (0)U (41) 


for an exponential model. 

D. Oblique Paths 

For an oblique path with elevation angle i9, we should replace Id with H esc i9. Amplitude variances 
for three profile models with a point receiver become 

(X 2 ) = 0.307C 2 fc 7/6 iI 11/6 (csc i?) 11/6 (Np) 2 for a slab model (42) 

( X 2 ) = 0.563C 2 fc 7/6 AffF 5 / 6 (csctf) 11/6 (Np) 2 for a thin layer model (43) 

( x 2 ) = 0.530C 2 o fc 7/6 ff 11/6 (csci?) 11/6 (Np) 2 for an exponential model (44) 
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V. Antenna Aperture-Averaging Effects 

For a receiving antenna with effective antenna radius a r , amplitude variance is an integration over the 
entire antenna area [2, 4, 7, 8]: 


<X 2 > = 


1 

A ? 




d 2 a 2 ix(Pi)x(P2)) 


1 




#1 


X 


(•2tt 


p 2 dp 2 / d<t > 2 {x{pu4>\)x{P2i fo)} 


(45) 


Figure 4 shows two surface elements, <J\ and a 2 , with a scalar distance |<Ti — 02 I on the receiving 
aperture. Spatial cross-correlation between two surface elements will generate a joint term Jq(k \ G\ — a 2 \) 
for the two-dimensional integration below. 


Based on the spatial covariance for a plane wave shown in Appendix B, the joint amplitude variances 
from two adjacent elements within an antenna dish can be expressed as 


MdM^)) = 47 J KdK$ n (K,)J 0 ( k \<Ti - a 2 1) J sin 2 dz 


poo pH 

27 r 2 fc 2 / ndnQ n (K) Jo ( k \(Ti — <r 2 \) 

Jo Jo 


1 — cos 


ZK 


dz 


poo 

27 r 2 Hk 2 / ndK& n (n) J 0 (k |<Ti — <r 2 1 ) 
Jo 


1 - 


sin(iiK 2 /A;) 
Hk 2 /k 


(46) 



Fig. 4. A two-dimensional integration is performed over a circular 
dish antenna with finite radius a r . Two surface elements show 
how they are cross-correlated, and the fluctuations received by 
surface elements are spatially averaged out due to the reduced 
gain factor for a large aperture antenna. 
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The spatial average over the entire antenna surface is 


<X 2 > = 


2 Hk 2 
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J JO Jo Jo 


(47) 


The addition theorem for Bessel functions separates the integrations as shown in Eq. (A-3). The 
angular integrations result in only one term ( n = 0) in the series left. Using the expression shown in 
Eq. (A-4), we have 


pd r p2tt pd r p 2'K / j ' 

J pidpi J dfii J p 2 dp 2 J dcj) 2 J 0 ( k\J p\ + p\- 2pipi cos(0 2 - &i) 



2Ji(na r ) 

)- 
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(48) 


Thus, 


(X 2 ) = 4tt 2 k 2 J ndnj dz$ n (K, z) sin 2 


2 Ji(/^a r ) 


K,a r 


(49) 


or 


poo pH 

(X 2 ) = 2ir 2 k 2 / ndn / dz$> n (n,z) 
Jo Jo 


1 ( zk2 


2 Ji(na r ) 


t\iCLr 


(50) 


where the term (2Ji(na r ) / na r ) 2 is the aperture-averaging wavenumber weighting function (or the Airy 
function) . It acts like a lowpass filter and eliminates the contribution from those turbulence eddies smaller 
than the radius of the receiver, a r . 

A. Slab Model 

Applying the Kolmogorov spectrum with a slab profile model of thickness H, we have 


(x 2 ) = 2n-Hk 2 / Kdn0.033C 2 n 


2 (n/ 3 ) ( _ sin {Hn 2 /k) 


Hn 2 /k 


2Ji(na r ) 
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= 0.651 k 2 HC 2 


dn 


Jo «' 


8/3 


1 - 
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2Ji(na r ) 
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= 0.307C^ 11/6 fc 7/6 Gi 



(Np) 2 


(51) 
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where a r ^/2Tr / HX is the ratio of antenna size, a r , over the first Fresnel zone size, \J H A/2-7T. After we use 
substitutions of 77 = a r \j2-K jHX = a r \Jk/H , £ = Hu 2 /k, the gain factor can be defined as 


Grfa) = 1.060 ^°°^6 



2Ji (WQ 
vVC 


(52) 


The integral G\{rj) cannot be solved analytically unless some approximations or simplifications are 
made. However, we can solve it by the numerical integration. The numerical solution for the gain factor 
is shown in Fig. 5 as a function of the ratio (rj) of antenna size over the first Fresnel zone size. If we 
normalize Eq. (51) by the point receiver solution, Eq. (33), we will have 


(x 2 M) 

<* 2 (o)> 


G 1(77) 


(53) 


The ratio of amplitude fluctuations between a finite aperture receiver and a point receiver is its gain 
factor. As expected, the antenna aperture plays a role in smoothing the amplitude fluctuations, because 
the gain factor, G{rj ), is always less than 1 and also rapidly decreases with increasing aperture size. Thus, 
the amplitude variations for a large antenna are always less than those for a point antenna. 



0.1 1 10 
r| = a r (k/H) 0 - 5 


Fig. 5. Gain factors for three types of turbulent layer 
models as a function of antenna size. The ITU model 
is also shown with a cutoff at G (rj) = 0.1. 
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B. Thin Layer 

Assuming a layer of thickness AH at altitude H and also taking AH << H, we have 


where 


<X 2 > 


0.563ClAHH 5/6 k 7/6 G 2 



(Np ) 2 


G 2 (r?) = 0.578 



^TJei l ~ C0S ^ 


2Ji(WC) 
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The gain factor for a thin layer model is plotted in Fig. 5. 

C. Exponential Model 

Similarly, we have 


(54) 


(55) 


where 


(X 2 ) 


0.530G 2 o ff 11/6 /c 7/6 G 3 



(Np ) 2 


G 3 ( V ) = 0.615 


f°° d( C 2 

'2Ji(rn/C)' 

Jo C 11/6 1 + C 2 

1 

**> 

sr* 

i 


(56) 


(57) 


The gain factor for an exponential model is also shown in Fig. 5. As compared with the other two models, 
the exponential model has the largest effect. We also see that all gain factors go to the unit when rj — > 0. 
The amplitude variance becomes the solution of a point receiver. The gain factor does not depend on 
either antenna size or Fresnel length alone. It strongly depends on the ratio (??) of antenna size over the 
first Fresnel length. 

D. Comparison with the International Telecommunication Union Model 

Crane and Blood 3 generate a piecewise function to express the antenna aperture-averaging gain factor 
for amplitude scintillation. This function defines the International Telecommunication Union (ITU) model 
as 


G{a r ) 


1.0 -1.4 (a r /V7TX^j 
< 0.5 — 0.4 
. 0.1 


for 0 < a r /VHX < 0.5 
for 0.5 < a r /VHX < 1.0 
for 1.0 < a r /V H\ 


(58) 


The results for this model are also shown in Fig. 5. We can see that the ITU model fits a slab model very 
well when a r /yf H\ < 1.4. 


3 R. K. Crane and D. W. Blood, Handbook for the Estimation of Microwave Propagation Effects-Link Calculations for 
Earth-Space Paths, Tech. Rep. 1, Document P.7376-TR1 (internal document), prepared for NASA Goddard Space Flight 
Center, Greenbelt, Maryland, June 1979. 
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VI. Scaling Factors 

Based on the above analytical solutions, we can obtain the following amplitude scintillation depen- 
dences on the wavelength, antenna size, and elevation angles. 

A. Frequency Dependence 

Noting the electromagnetic wavenumber is k = 27r/A, we have for a point receiver 


„ (K \ 7/6 

(xl) V fc 2/ \Xl) V/2/ 


(59) 


when the signal’s wavelength and frequency have the relation of A / = v g and assuming that the wave 
propagation speed, v g , is a constant. 

For example, relative to X-band (8.4 GHz), the amplitude variance at Ka-band (32 GHz) will increase 
by a factor of 4.8. 

The gain factor for a finite aperture antenna expressed as a function of wavelength is 


B. Antenna Size 

Normalizing Eqs. (51) through (53) by Eqs. (33), (36), and (41), respectively, we have 

(x 2 K)) _ r ( / ~2tT 

(X 2 (0)) ^ r V H\ 


(60) 


(61) 


which also depends on what type of turbulent profile models one chooses. 

C. Elevation Angle 

From Expressions (42) through (44), we can find 


(x 2 C#i)) 

<XW> 



11/6 

for > 5 deg 


(62) 


The relation will be broken up at 5 deg and below because of multipath effects at low elevation angles. 


VII. DSN Goldstone Receiver 

Now we can apply these analytic solutions to a 34-m parabolic dish antenna at the Goldstone DSN 
site at Ka-band. We will use the following parameters for the Goldstone site calculation [9,10]: 
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C 2 n0 = 0.5 x 10- 13 (m-( 2 / 3 )) 


A = 0.01 (m) 

H = 8.0 x 10 3 (m) 
a r = 34/2 x 55% (m) 


Here we have assumed that the dish antenna has a 55 percent efficiency relative to its physical radius. We 
will show an example of how to calculate the amplitude variations, as follows. From Fig. 6, we can find 
that, when 77 = a r ^2n/XH = 2.62, G\(if) = 0.075 for a slab model; = 0.12 for a thin layer model; 

and Gs(rf) = 0.19 for an exponential model. Thus, when elevation angle •& = 90 deg, for an exponential 
model we have, from Eq. (56), 


(x 2 ) = 0 .530 C 2 o 11/6 /c 7 / 6 (sin 1?) " ( n / 6 ) G 3 



(Np) 2 


(63) 


<X 2 ) = 1.33 x IQ" 4 (Np) 2 


(64) 


Its root of mean square (rms), that is, amplitude fluctuation, is 

Xrms = 1.15 X IQ" 2 (Np) 


Using Eq. (6) and converting to decibels, we have 

Xrms = 0.05 (dB) 
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Fig. 6. Elevation-angle dependence of amplitude fluctuation (in 
decibels) for both a 34-m DSN antenna and a point receiver for 
three types of atmospheric profile models. 
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When elevation angle $ = 20 deg, 


<X 2 ) = 9.51 x 1CT 4 (Np) 2 
Xrms = 3.08 X 1(T 2 (Np) 


Xrms = 0.134 (dB) 


Figure 6 shows the amplitude fluctuations for elevation angles from 20 deg to 90 deg for three atmo- 
spheric turbulent models for a DSN 34-m antenna at Goldstone; the DSN antenna will not operate below 
the 20-deg elevation angle for Ka-band. For comparison, the amplitude variations for a point receiver 
also are shown in the figure. To calculate the amplitude fluctuation for a thin layer, we have assumed a 
400-m layer thickness. Of the three models, the exponential model has high amplitude variations because 
of its large gain factor. The ratio of amplitude changes between a receiver with finite radius a r and a 
point receiver is 


Xrms xP'r , 
Xrms (0) 



(65) 


VIII. Summary 

On the basis of the Kolmogorov turbulence theory, we have derived all analytical solutions for the 
amplitude variances for a finite aperture antenna for the DSN scenario. Through this analytical study, we 
found that amplitude scintillation is an important factor that cannot be neglected under the clear weather 
at Ka-band, even for the large antenna used at Goldstone. At low elevation angles, the scintillation 
intensity remains high. For example, at a 20-deg elevation, amplitude fluctuation can be as much as 
0.13 dB for a 34-m DSN antenna for an exponential turbulent model. If these fast fadings are further 
superimposed on some slow fadings (caused by water vapor density changes), deep fading will definitely 
cause a problem for the DSN link. Using these fluctuation values, we can further study power spectra 
and fading rates of the amplitude scintillations [11]. 

A receiver antenna with finite aperture will have different effects from a point receiver on the scin- 
tillation intensity and temporal characteristics. Because the receiver’s output is a spatial average of the 
wavefront fluctuations over the aperture, this results in the observed fluctuations being smaller than 
those received by a point receiver. We can see that the gain factor for a finite aperture antenna strongly 
depends on the ratio (rf) of antenna size over the first Fresnel zone size. The gain factor also is dependent 
on what types of atmospheric turbulence profile models we use. 

These results should help in system design and fading prediction for Ka-band downlink for a low- 
margin telecommunication system such as the DSN. It also establishes the theoretical basis for any future 
experimental comparison at Ka-band. The fading study on atmospheric scintillation also provides us a 
possible mitigation approach in controlling the uplink power based on predicted amplitude fluctuations. 
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Appendix A 

Bessel Functions and Other integral Equations 


The Bessel functions are as follows: 


/ d(j>ex p [ rn r cos (<t> — w)] = 2nJ 0 (n r r) 

Jo 

f°° 2 1 f a 2 \ 

/ Jn(ax) cos(bx 2 )xdx = — sin ( — I 

Jo 2 b \4b) 

J 0 (k\J p\ + p 2 ~ 2 ~pipi COS (02 - 4>l) ^ = T. S n Jn(Kpi)J n (np 2 ) COS [n(<j > 2 ~ 4>l)] 


i r 

— ^ / xJo(xn)dx = 

a 2 Jo 


Ji(na) 


na 


The Delta functions are as follows: 


pR pR pR 

/ dx dx'f(x)g(x') exp [m x (x — a/)] =27 tS(r x ) / dxf(x)g(x) 
lo Jo Jo 


5(rcos'i/j) = —5 (cos 'll;) 

R 


dtp f (VO 5 (k cos ip) = -f C 11 

K \ 


K \2 


The integral solutions are as follows: 


r <k 
Jo C 11/6 



0.943 




[1 -cosC] = 1-729 



x m l dx 
1 + x n 



(A-l) 

(A-2) 

(A-3) 

(A-4) 

(A-5) 

(A-6) 

(A-7) 

(A-8) 
(A-9) 
(A- 10) 
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Appendix B 

Antenna Aperture-Averaging Equation 


The amplitude fluctuations for a plane wave measured at separated locations can be described by 

X(R) = —2 k 2 J d 3 rAi y R,r)5n{r) (B 

X(R + A p) = — 2 k 2 J d 3 rA(R + A p, r)5n{r ) 

= - 2 k 2 J d 3 r'A(R, r')6n(r' + A p) (B 

where A p = p 2 — Pi is a vector in a plane of a receiver antenna: 

(x(R)x(R + A />)) = 4 k A J d 3 rA(R, r) J d 3 r'A(R , r) ( Sn(r)6n(r 1 + A p)) (B 


where 


( 5n(r , t)Sn(r' + A p, t)) = J d 3 K$ n (n) exp [in ■ (r — r' — Ap)] 


(x{R)x{R + Ap)) = 4/c 4 J d 3 i'A(R } r) exp(i/c • r) J d 3 r'A(R, r’) exp(*/c • r') 


x J d 3 K^ n (n) exp (in • Ap) 


The three-dimensional integration on k is written in spherical coordinates as 


nO O n OO nix nZTT 

/ d/-c 3 exp(i/c • Ap) = / K 2 dn dip sin ip5(n cos ip) / cko exp(in ■ Ap) 

Jo Jo Jo Jo 


In the usual case, ip — tt/2 and sin^ = 1, so that 


nOO n OO n2lV n OO 

/ cIk 3 exp('</c • Ap) = / / dwexp(*/cApcosw) = 27 t / ndnJo^Ap) 

Jo Jo Jo Jo 


because 
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4 k 4 J d 3 rA(R,r) exp(i/c ■ r) J d 3 r'A(R, r') exp(i/c • r') = 47r 2 /c 2 J dzsin^^—^j 


(x(R)x(R + Ap)) = 4n 2 k 2 J dz sin J dKK$ n {n) J 0 {nAp) (B-6) 

An additional term Jo(nAp) that spatially joins together two elements separated by A p is included in 
the spatial covariance. For the antenna aperture averaging, we have 

(xOi)x(> 2 )) = 47r 2 fc 2 y ndn$ n {K)J 0 {K |cti - a 2 \) J sm 2 dz (B-7) 


where a \ and 02 are, respectively, two elements within the antenna. 
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